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This article generalizes the notion of the local density of a many-body system to introduce col- 
lective coordinates as explicit degrees of freedom. It is shown that the energy of the system can 
be expressed as a functional of this object. The latter can in turn be factorized as the product of 
the square of a collective wave function and a normalized collective-coordinate-dependent density. 
Energy minimization translates into a set of coupled equations, i.e. a local Schrodinger equation 
for the collective wave function and a set of Kohn-Sham equations for optimizing the normalized 
density at each point in the collective space. These equations reformulate the many-body problem 
exactly provided one is able to determine density- and collective-wave-function-dependent terms 
of the collective mass and potential which play a similar role to the exchange-correlation term in 
electronic Kohn-Sham density functional theory. 

PACS numbers: 21.60.De, 21.60.Jz 
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I. INTRODUCTION 

Dcnsity-functional-based models, being the only mi- 
croscopic, fully quantum-mechanical tool currently avail- 
able to provide insight on the structure of nuclei up to 
the heaviest ones, are the focus of intense investigation to 
improve their accuracy and precision. For example, re- 
cent developments in effective field theory methods and 
the theoretical foundations of nuclear DFT have triggered 
attempts to bridge ab-initio many-body methods, based 
on chiral effective field theory and high-precision two- 
and three-nucleon interactions, and DFT, using the for- 
mer to build parts of the latter Another current 
line of work consists in enriching the form of empirical 
energy functionals in a systematic way (9j-[Tl|. 

One puzzling fact about nuclear DFT is that it owes 
much of its power to the use of symmetry-breaking den- 
sity configurations. The method thus deviates from the 
symmetry-conserving Hohenberg-Kohn (HK) [l2j and 
Kohn-Sham (KS) [13j framework. The most basic exam- 
ple of a broken symmetry is the translational invariance 
of the Hamiltonian, which is troublesome for self-bound 
finite systems best described by a localized density. KS 
frameworks have recently been built for a trapped sys- 
tem 1J] and for a functional of the internal density, i.e., 



expressed in the reference frame of the nuclear center-of- 
mass EE [l~o| . A density- functional framework allowing 
to break arbitrary symmetries was put forward in Rcf. 
[l7j , relying however on an unspecified restriction of the 
variational space for a trial wave function. Another ex- 
ample of a broken symmetry is the non-conservation of 
particle number in the BCS treatment of pairing. As an 
alternative, exact solutions of the pairing Hamiltonian 



and symmetry-restored Hartree-Fock-Bogolyubov equa- 
tions have been formulated as a functional (function, in 
fact) of occupation numbers in a given single-particle ba- 
sis [13t20|: and functions of occupation numbers have 
been studied as an alternate route to building mass tables 
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21| . These works keep the HK/KS formalism unchanged 
except for a variation in its basic degree of freedom: the 
system's single-particle density. 

However, formally integrating into DFT the break- 
ing of rotational invariance, as well as the violation of 
particle-number conservation, remains to be done. Al- 
though effort is currently being put into building func- 
tionals of the scalar, symmetry-conserving density [22J, 
it is likely that breaking these symmetries is essential to 
describe what is understood as nuclear deformation, and 
treat pairing in finite systems, while keeping computa- 
tional complexity to a minimum. 

The cost-efficiency of KS-like schemes comes from the 
idea of using a Slater determinant to reproduce the lo- 
cal density p of a given correlated wave function, thereby 
capturing essential quantum effects, and encode the miss- 
ing correlations into a functional E xc [p] @,|2iJ. Thus, no 
explicit mention is made of a many-body wave function. 
Nonetheless, current expectations about nuclear DFT 
were mostly raised by the success encountered during its 
first life as mean-field theory performed with effective 
density-dependent interactions, together with beyond- 
mean-field extensions such as the Generator Coordinate 
Method (GCM) performed on top of symmetry-restored 
mean-field states [2oT - |28j . This was commonly under- 
stood as an approximate wave function method, which 
requires the restoration of broken symmetries for con- 
sistency with the underlying Hamiltonian, and allows 
multiple symmetry-restored configurations to be mixed 
to describe zero-point collective motion, or shape coex- 
istence, and make extensive spectroscopic predictions in 
nuclei where these degrees of freedom are important. The 
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somewhat ad-hoc extension from a Hamiltonian picture 
to a dcnsity-functional-based one leads to pathologies in 
the theory [29[ , which can be formally addressed (20, Hl| 
at the expense of additional complexity and constraints 
on the form of the functional, such as forbidding terms 
other than polynomials of the density [32l - [34l |. This is 
referred to as the Multi-Reference Energy Density Func- 
tional (MR-EDF) model. 

An alternate formalism uses a collective "Bohr" Hamil- 
tonian (351442^ , initially developed as a model for a quan- 
tum vibrating liquid droplet, as an alternative or approx- 
imation to the Hill- Wheeler equations of the GCM. The 
mean-field energy landscape in the space of deformation 
parameters is then used as a potential (with or with- 
out zero-point energy corrections) and mass parameters 
determined from the Slater determinants enter a corre- 
sponding kinetic operator. 

MR-EDF and collective Hamiltonians are powerful 
tools, well adapted to the idiosyncrasies of the nuclear 
many-body problem, yet they do not have a clear con- 
nection to the first-principles formulation thereof, which 
may limit their future development. Here, I attempt to 
improve on this situation. 

The present article aims to formulate a density func- 
tional theory allowing to break spatial symmetries (e.g. 
translation, rotation) and treat collective motion in 
a spirit similar to the GCM, MR-EDF, or collective 
Hamiltonian-bascd methods, while keeping the theory ex- 
act in the DFT sense, i.e. provided we can determine 
the exact functional. First, the relevant mathematical 
objects are defined, then the existence of a functional 
is summarily proven in Section [II] A useful form of the 
latter is given in Section IIIII The example of transla- 
tional motion is used as an illustration in Section IIVI 
Finally, a minimal Kohn-Sham-like scheme for introduc- 
ing auxiliary Slater determinants is presented in Section 
IVl Results are discussed throughout and summarized in 
Section EH 



II. DEFINITIONS AND EXISTENCE OF A 
FUNCTIONAL 

Consider a system of N fermions with Hamiltonian 
H = T + U + V, where T, U and V are respectively 
kinetic, interaction and external potential terms. Trial 
A-body antisymmetric wave functions ^ depend on N 
coordinate 3-vectors fi , i = 1 . . . N (omitting, in this first 
presentation, internal degrees of freedom of the particles 
for the sake of clarity). Let us write 



T=-±A = (3) 

i 

f(R)=]>> cx t(r 4 ), (4) 

i 

while, for now, we shall keep U as an arbitrary inter- 
action term. 

Now consider a set of diffcrcntiable real functions 
Qn(r) indexed by /i = 1 . . .n. We will use these as po- 
tential operators, 

Q^R) = Y,Q^), (5) 

i 

their expectation values q^ in a many-body wave func- 
tion defining a set of collective coordinates. Now, for any 
given set of values q = (gi, . . . , q n ), each being taken 
in the interval of possible values of (R) , we can define 
the operator 

P(q,R) = J] ^(R) (6) 

This operator selects configurations of the N particles 
for which the collective coordinates defined by the func- 
tions Q/x correspond exactly to the given values (q). It 

effectively projects ^ onto an eigenspace of the Q^. Using 
the definition above, it is trivial to prove the projector- 
like property 

P(q, R) P( g ', R) = 5^(q - J) P(q, R). (7) 
as well as the closure relation 

j d n qP(q,R) = l, (8) 

where the g^-integral runs, as in the following, over the 
interval of possible values of Q P (R). We can now use P 
to define the generalized density, 

D(q, r) = N J d 3N R c5 (3) ^- n) P{q, R) 

x tf*(R) \P(R). (9) 
Comparing with the usual density, 



R=(f 1 ,...,r JV ), (1) 
d 3Ar REEd 3 r 1 ...dV Jv . (2) 

The kinetic and external-potential terms of H can be 
written as (in units where h = m = 1) 



p(f f ) = N J d 3N RS( 3 \r~r 1 ) **(R) *(R), (10) 

we see that D(q, f) is, up to a factor N, the proba- 
bility density of finding one particle at r and the collec- 
tive configuration of the N particles at q. As P is an 



3 



iV-body operator, D contains information from up to N- 
body components of the density matrix associated with 

It is possible to split this information by defining, first, 
a collective wave function (cwf) 



f(q) = e 



16(g) 



1 J d 3 rD{q,r) 



1/2 



(11) 



rf 3JV R$*(R)P(g,R)$(R; 



-,1/2 



,(12) 



where 6(q) is a chosen phase depending on the prob- 
lem at hand. We shall assume there exists a natural and 
unambiguous choice. For ground states, Q(q) = seems 
appropriate; otherwise, an irreducible representation of 
a symmetry group of H in g-space may provide the de- 
pendence of 9(q) on some or all coordinates (this will be 
further discussed Sections Mil and HV|) . 

Second, we can define the g-dependent density (defined 
first almost everywhere, then elsewhere by continuity), 



d{q,r) = \f{q)\- 2 D{q,f), 



(13) 



which captures the conditional probability density of 
finding a particle at r if the collective configuration is 
q. The meaning of d(q,r) can be made more explicit by 
introducing the "slice" wave function 

v%R) = /-%)P(g,R)vf.(R), (14) 
which satisfies, using Eqs. ([7|) and (fT2j) 



J d 3Ar R<T(g,R)*(g',R) = 8 {n \q-j), (15) 

fd n qf(q)*(q,R) = *(R). (16) 

This wave function is non-zero on a manifold of dimen- 
sion 3iV — n determined by the values of the collective co- 
ordinates. Then, d(q, r) is the density of this state (using 
Eq. again), 

5 { - n) {q-j)d{q,r)=N J d 3N R 6 {3) {f - n) 

x <T(g,R)*(g',R). (17) 

Moreover, the quantities defined above satisfy the nor- 
malization relations, obtained from their respective defi- 
nitions and the closure relation, Eq. (|5J), 



N, 


(18) 


1. 


(19) 


N, 


(20) 



and the following relations, 

J d n qD(q,f)=p(r), 
d 3 rD( qi r) = N\f(q)\ 2 , 



(21) 
(22) 



which exhibit the role of D as a joint probability dis- 
tribution and those of p and |/| 2 as the corresponding 
marginal distributions. Finally, we can see that the value 
of q is encoded in d(q, r) through 



d 3 fQ )i (r)d{q,r) 



(23) 



\f(g)r J d R Q M (R) P (g , R) 

x**(R)*(R) = g M , (24) 

thus reducing the variational domain for d. We could 
have used a two-body or higher operator for Q, but this 
would prevent us from obtaining a simple expression like 
Eq. ([23]) for this purpose. 

Now, let us see how the energy of the system can be 
expressed as a functional of D. Let w(q, r) be a real 
function bounded from below and 

w(q,R)=J2 w (g^i), (25) 

i 

W(R) = [ d n qw(q,R)P(q,R). (26) 



This N-body operator applies on each slice of a wave 
function a different, (/-dependent, local, single-particle 
potential w(q, r): 

W(R) *(R) = J d n q /(g) w(q, R) R). (27) 
Its expectation value is, using Eqs. (fTB"]) , (fT7|) and (|27p. 



|^| = J d n qd n q; ridfW) 

x J d 3N R R')w(g',R)*(g',R),(28) 
d n q I d 3 fw(q,r)D{q,r), (29) 



i.e. this potential is purely multiplicative with respect 
to the generalized density D, just as a local one-body 
potential v(r) is with respect to p(r). In fact, such a one- 
body potential is a special case of W with no dependence 
on q (w(q,r) = v(r) in Eq. ([25])). 

Let us add such a potential to the Hamiltonian and 
define the following functional, which implies solving for 
the ground state of H + W: 
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F[w] 



(30) 



To obtain a physically useful theory, we first need to 
ensure that w and ^, hence all physical observables, are 
functionals of D. The proof is identical to the usual case 
[HI EH Hi| and shall not be repeated here. 

Using the Hcllmann-Fcynman theorem and Eq. (|2T)|) . 
we have 



6F[w] 
Sw(q, r) 



D(q,r). 



(31) 



We can thus use a Lcgcndre transform to write a func- 
tional of D, 



E[D] = min 



F[w] - J d n q J d 3 rw(q,r)D(q,r) 



= min (^\H\ 1 <S>) 



(32) 
(33) 



where ^ — > D means that the variational domain for 
*5f is restricted to wave functions having the generalized 
density D. It thus appears that the energy of the system 
can be written as a functional of D, in a similar fashion 
to the HK result. The ground-state energy of the system 
can then be found by minimizing the functional E[D\. 

Since D depends on more variables than just the coor- 
dinates of one particle, it introduces additional degrees of 
freedom compared to standard DFT. Moreover, as seen in 
SectionHH given a phase choice, D can be unambiguously 
decomposed into a CWF / and a g-dependent normalized 
density d. We can thus also write 



E[D] = E[f,d\. 



(34) 



This alternate formulation will become useful when de- 
riving the formal basis for a practical many-body method 
in the next few sections. For now, let me give a few 
examples of collective coordinates that can be usefully 
incorporated. 

For instance, assuming a vanishing external potential 
(V = 0), using x/N, y/N and z/N for Q^{r) (x,y,z 
being the Cartesian components of r), the collective co- 
ordinates q^ are the components of the center-of-mass 
(CoM) coordinate vector R. The generalized density 
is then D(R,r) = \f(R)\ 2 d(R,r), where, provided U is 
translation-invariant, / = O -1 / 2 with Q a normalization 
volume, d(R,r) = /Oi nt (r — R), and p; n t is the internal 
density of the system, non- vanishing if the system is self- 
bound. In this case, pint is the only physically relevant 
degree of freedom present in D; the energy is a functional 
of the internal density, which is the result of Ref. [l6| that 
we recover here as a limit case. 

More generally, when the cwf f{q) is known from sym- 
metry arguments and the coordinate-dependent density 



for one value of the coordinates, d(q, r), can be deduced 
from its value at some natural reference point d(0, r) by 
a symmetry transformation, the latter is obviously a suf- 
ficient degree of freedom. 



Suppose we now add the functions (y + z 



x 2 + z 2 , 



x +y , —xy, —yz, —zx) to the set of Q M (r ) . We can now 
use as coordinates the components of the inertia tensor 
of the nucleus in the laboratory frame, 



<Tr 



y 2 + z 2 
-yx 
—zx 



-xy 
+ z- 
-zy 



x 2 + z 2 



-yz 



d(RJ,', 



(35) 



From this and the CoM position, using the Huygens- 
Steiner theorem from solid mechanics, we can recover the 
inertia tensor in the CoM frame Jo, 



fa = J-N(I?I-R®R), 



(36) 



where I is the identity matrix and ® is the tensor 
(outer) product of vectors. In turn, Jo can be trans- 
lated into a root-mean-square matter radius (r rms = 
Tr(J )/2iV) and a quintuplet of quadrupole deformation 
parameters and Euler angles (/?, 7, <p, ip) [HI, [4l|. We 
can thus treat deformation and rotation degrees of free- 
dom dynamically. 

At this point we might find that the number of col- 
lective coordinates in our theory, i.e. eight independent 
components of (R, Jo) is actually too large for practical 
applications. For example, fluctuations of the radius are 
not usually considered an essential dynamical degree of 
freedom. We can define a functional with coordinates 
removed as follows. Consider a reduced set q of n < n 

coordinates, and the generalized density D(q,r) depend- 
ing on this reduced set. The energy can be written as 



E[D] = min E[D], 

D^D 

where D — > D means that for all r. 



dqn+i ...dq n D(q, r) = D(q, r). 



(37) 



(38) 



This operation can be performed after a nonlinear 
transformation among the coordinates q, such as the one 
mentioned above to obtain the canonical Bohr coordi- 
nates from the inertia tensor. Using this, we can remove 
in succession r rms , then, if desired, 7 and ip (to obtain 
a functional describing only axially-symmetric deforma- 
tion). Alternatively, we may remove (3 and 7 and keep 
only Euler angles as coordinates. In this case, again, 
the cwf as well as the transformations of the density 
d(R,ip, , d,ip;r) are known analytically, and we can ex- 
press the energy as a function of a single, deformed in- 
trinsic density. 
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Note that Eq. (J55J) has a similar form to Eq. (l2Tj) . 
We could, in principle, use this procedure starting with 
the full TV-body local density **(R)*(R) and, inte- 
grating out collective coordinates, yield a succession of 
generalized-density functionals, all the way down to the 
Hohenberg-Kohn functional of p(r) if we remove all of 
them. This formalism thus appears very general and flex- 
ible. In particular, it gives us a choice between a "single- 
reference" description of many-body systems in terms of 
a single density (if symmetries define the dependence on 
q of the relevant quantities), of a "multi-reference" de- 
scription explicitly coupling single-particle and collective 
motion. 

Another important point to stress is that this theory 
is symmetry-conserving. Let S(r) be a an orthogonal 
transformation of the coordinates. Then, define 

5(R) - (S(n),S(r 2 ),...,S(r N )), (39) 

and suppose the set of collective coordinates is chosen 
so that we can define 

S(Q(R)) = Q(S(R)). (40) 

Then, the projector P(q, R) is invariant under the si- 
multaneous transformation of single-particle and collec- 
tive coordinates, 

P (S(q), 5(R)) = 6™ (Q(«S(R)) - S(g)) , (41) 
= <5W(Q(R)-g). (42) 

If, furthermore, the transformation leaves the many- 
body wave function invariant up to a phase 77, 



detailed discussion, which involves the phase choice en- 
tering the definition of / and needs to be done on a case- 
by-case basis. 

One remaining hurdle, is pairing: treating superfluid 
systems by breaking the conservation of particle number 
to yield a non-zero pair density would require an opera- 
tor P(q) projecting on the associated U(l) gauge angle 
(yielding a particle-number breaking slice \^>(q))). This, 
in turn, requires the definition of a gauge-angle operator; 
technicalities of such phase operators have been worked 
out in the field of quantum optics |45| starting from a 
particle-number representation. In our case, this would 
entail a switch to a Fock-space representation, which is 
cnvisionable yet cumbersome (the present derivation is 
pervaded with iV-body operators), and, moreover, the 
definition of a particle-number basis corresponding to 
each point (or state) in the iV-body Hilbert space, i.e. 
choosing one definite particle-addition operator. Such 
a procedure relies on an arbitrary choice (what is the 
N + 2-body component of the iV-particle system's corre- 
lated wave function ?) that has to be studied in detail, 
and is beyond the scope of the present work, the remain- 
der of which shall be concerned with normal systems. 
Let me simply suggest that pairing should probably be 
treated with a different, simpler scheme, which shall be 
described in a future paper. 



III. COLLECTIVE SCHRODINGER EQUATION 

The rest of this article shall be devoted to deriving the 
formal basis for a practical many-body method based on 
the result from the previous section. Here I focus on the 
dependence of the energy on the cwf / and derive a form 
of the functional that allows to optimize the latter. Let 
us define a trial wave function 



*(S(R))=»?*(R), (43) 
we have, from the definition of D, 

£>(5( £ ),5(f)) =D{g,r f ). (44) 

This simultaneous transformation thus leaves D invari- 
ant as well, the same property being obtained for d if / is 
invariant under S_. Of course, for a given, fixed value of q, 
D{q 1 r) and d(q, r), understood as functions of r alone, do 
not have to be symmetry- invariant. This brings a solu- 
tion to the usual conundrum around using the symmetry- 
conserving HK framework to justify nuclear DFT. 

For completeness, the case of a wave function trans- 
forming as a non-trivial representation of a symmetry 
group should be mentionned. In general, D is not in- 
variant in this case. Ideally, the collective coordinates 
should be chosen so as to replicate the group structure 
in the transformations of / and d. This deserves a more 



V = J d n gg(q)\y(q)), (45) 

where g is a trial cwf, whereas in the following /, d 
and #(£, R) are defined through Eqs. (O, Jl3|) and 
from a starting wave function \1/(R). This will allow us 
to examine the dependence of the energy on the cwf and 
put E[g, d] in a convenient form, then set g — f at the 
end. The trial energy is 

E=(^'\f + U + V\^'), (46) 
= J <Tqd n i g*{$g{£){*{£\f + U + V\*{£)), 

(47) 

Let us first give the interaction matrix element. Here 
we assume a local interaction (which can nevertheless 
contain three-body or higher operators), 
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(tt(g)|#|*(g , )> = 



(P(g)) 

xP(g,R) i7(R)tf*(R)tf(R), (48) 



with a similar expression for V', where we define 

(P(q)) =J rf 3iv R%,R)r(R)f(R) (49) 
= l/(9)| 2 - (50) 



The kinetic matrix element, in turn, is (after integrat- 
ing by parts; here and below we assume Dirichlet bound- 
ary conditions or an infinite integration domain, allowing 
to drop boundary terms) 



<*(g)|T |v%')> = -J d 3N R 

x Vtf*(g,R) • V#(g',R), (51) 



d 3N H 



2 /*(£)/(£') 
x{v%R).VP(g',R) 
x **(R)*(R) 
+ P(g,R) VP(£',R) 

• V**(R)*(R) 

- **(H) Vtf(R) 

-5(")( £ -g')P(£,R) 
x tf*(R) Atf(R)}. 



(52) 



The gradient of P can be derived using Eq. J6]) and the 
chain rule, 



VP( £ , R) = - J] VQ M (R) a M P(£, R), (53) 



/if 

xM/ _1 (g)s(9)) 
-*E'V*~W(£) 

xd ll {f- 1 (g)g(q)) 
- J d m R P{q, R) **(R) A *(R) 

xr^ 1 (£).9*(£)r 1 (£).9(£)},(54) 

where we introduce 



(55) 



(56) 



f ,v RP(g, R) 

x VQ M (R) • VQ^R) **(R) *(R), 



rRPf?, R)VQ M (R) 
V$*(R) *(R) - **(R) V*(R) 



Using Eq. (|53[) and the steady-state continuity equa- 
tion for the probability current of ^(R), we can easily 



check that J M itself satisfies 



(57) 



Its definition and this property suggests its role as a 
collective current. 

Finally, replacing / everywhere by its expression in- 
volving 9(q) and (P(q)), using integration by parts again, 
then setting g = /, we have 



E[f, d]=J d n q f* (g) ~Y1 d f A ^ (£) 3 " + U ^ 



/(g), 



(58) 



where we made act on all factors to its right. The 
potentials entering Eq. ([55)) are defined as follows. First, 
the collective mass term is 



where indicates differentiation with respect to q^. 
Applying the result of Eq. ([52]) in Eq. (|47|) . using inte- 
gration by parts to transfer on the cwf, and reducing 
the double integral with Eq. ([7]), we can write the kinetic 
energy as 



(>(«)>' 



(59) 



This collective-mass term does not explicitly involve 
the interaction; it consists of a part of the kinetic energy. 
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This is expected since a local potential does not couple 
slices with different q, the latter being non-zero on differ- 
ent, non-overlapping manifolds in the many-body coordi- 
nate space: a local operator only contributes to the local 
collective potential. The dependence of on the in- 
teraction is thus implicit, due to dependence on the wave 
function itself in (P(q)) and F )J , l ,(q). 

To the contrary, the GCM and its collcctive- 
Hamiltonian based approximations typically use, as 
building blocks, Slater determinants which are not lo- 
calized in collective coordinate space (hence their non- 
orthogonality and, in the GCM, the need to remove zero- 
norm states when solving the Hill- Wheeler equations). 
Our collective mass and the one found in collective- 
Hamiltonian models may thus have slightly different 
meanings. 

Second, the collective potential is 



U(q) = J2 F ^(g) 



ia M 0(gW(g) 

2 (P(q)) 



ld^(P(q))d„(P(q)) 



(P(qW 



+ z5> 



//// 



(P(q)) 



^d v (P(q)) 



d 3N RP(q,R) 



(P(g)) 

x **(R) [ - |A + U(R)] *(R) 
+ / d 3 rv xt{r)d(q,r). 



(60) 



If is an cigenstatc of H, the last two terms of Eq. ([60]) 
boil down to the energy E. In fact, all terms proportional 
to Ffa, and J M entering Eq. (|58|) cancel each other, as is 
made obvious by setting / = g prematurely in Eq. (1541) . 
and these expressions could be simplified to a trivial form. 
However this particular separation of the energy is useful 
in isolating the dynamics of the system with respect to 
the chosen coordinates, while integrating out uninterest- 
ing ones. We shall see below, with an example, that it 
yields a meaningful physical value for the collective mass 
and potential. 

This potential contains a piece of the kinetic energy, 
as well as all interaction and external-potential terms 
of the Hamiltonian (the latter being contained in the 
last line). Here, we assumed a local interaction term. 
A non-local one would simply make the potential itself 
non-local in the collective space, i.e. U(q,q'). Such 
a non-local interaction is commonly found as the re- 
sult of a renormalization-group (RG) evolution (46l - |48j 
of a starting, local model of the nucleon-nucleon inter- 
action. This suggests that our collective Hamiltonian 



is not renormalization-scale- invariant. Since D{q,r) in- 
volves components of many-body density matrices of the 
system, it is sensitive to details of the wave function and 
should not be considered an observable in the RG sense, 
or at best a scheme-dependent one , and the same has 
to be deduced for quantities entering Eq. ([55)). However, 
the main purpose of the present formalism is the descrip- 
tion of low-energy collective states and obscrvablcs which 
should not be sensitive to such details. The generalized 
density, as the cwf, should thus be largely scale-invariant 
in practice, for appropriate choices of the collective coor- 
dinates. We are thus presented with a scale-dependent 
collective Hamiltonian with largely scale-independent so- 
lutions, indicating that scale-dependence mainly occurs 
through reshuffling of contributions to the energy be- 
tween non-locality in U and the collective mass term. 
In practice, it should be safe to limit ourselves to param- 
eterizations of a local collective potential. 

Finally, the potential multiplying the current operator 



is 



d v 6(q) , Mq) 



(P(q)) (P(q)) 



(61) 



Here, the first term in the definition of V, as well as 
the term involving J M in Eq. (j60[) are proportional to the 
derivative of the phase introduced in Eq. ([lip and ensure 
invariance of the energy with respect to gauge transfor- 
mations of /. The notion that the energy (especially the 
kinetic energy) should not depend on the complex phase 
of the wave function may be counter-intuitive. It is worth 
reminding here that we are dealing with a functional of 
a generalized density, and that the cwf /, if useful for 
formulating a theory of collective motion, has been in- 
troduced somewhat artificially. Thinking in terms of the 
correlated many-body wave function of the system, the 
functional E[D] = E[f, d] yields the lowest energy of all 
states having D as their generalized density; the state 
which minimizes this energy is unique and we can ac- 
cess no other. For example, in the case of translational 
motion examined above, plane waves with non-zero mo- 
mentum are excluded from the theory, as they have the 
same density, in terms of the CoM coordinate vector, as 
the zero-momentum state. 

The second term in Eq. (p?l"[) involves the current J^{q). 
We can use the freedom of choosing the phase 9(q) men- 
tioned above to cancel it with the first term if we can 
make the phase satisfy 



J2 F ^(i) d » 9 k) = - J ^q), 



(62) 



which amounts to introducing in the cwf the current 
which is present in the underlying many-body wave func- 
tion. Cancelling this term is useful to simplify the func- 
tional into a form that is more convenient to later derive 
a Schrodingcr equation for /, 
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E[f,d\= / d n qf*(q) 



i a' 



fig)- 

(63) 



The phase will then simply come out from the solution 
to that equation. 

The collective inverse mass A fiv (q) and the collective 
potential U(q) depend on the wave function ty: they 
are, for each q, functionals of / and d. A collective 
Schrodinger equation can be obtained by minimizing 
E with respect to /, with a constraint on the norm, 
Eq. OH, viz. 



VQ„(R) ■ VQ„(R) = 5 -f, 

F ^(g) = Jf( p (g)) 



(66) 
(67) 
(68) 



The mass of the system entering the kinetic term is 
thus N times the mass of the constituent particle itself, 
as expected. 

Assuming w ex t(^0 = ana - ^ invariant under trans- 
lations and Galilei transformations, we can write the 
Hamiltonian as 



5{E[f,d]-E> 
Sf* 



= = 



_ o ^2 d ^u(g)d u 



+ U(q)+U Ia ,(g)-E'\f(q), 

where acts on all factors to its right, with a rear- 
rangement potential that appears because of the func- 
tional derivation of A and U, 



U la (q) ee / dV/V) 



<5|/| 2 (<?) 



(64) 



/(<?') (65) 



where 9' differentiates with respect to q'^ . Note that in 
general, a distinction has to be made between E (energy) 
and E' (eigenvalue of the collective Hamiltonian) due to 
the rearrangement energy. 

The remaining issue is then to optimise the q- 
dependent density d(q,r). This will be dealt with in 
Section |Vl 



(69) 
(70) 



and eigenstate wave function of the system as 



*(R) = * cm (i2)* tat (a), 



(71) 



where S is a vector of Jacobi coordinates allowing to 
describe internal motion of the particles. 

For D(R,r) independent of R, the value of ^ cm (R) 
that minimizes the cm. kinetic energy is ^ cm (R) = 
q-i/2_ r£ Q access a state with non- vanishing cm. kinetic 
energy, we need to choose a trial f(R) such that the cm.- 
coordinate density \f\ 2 (R) is inhomogeneous. One such 
choice is 



f(R) 



sm(K ■ R). 



(72) 



For d, let us set d{R,r) = p ln t(f), where p ln t is the 
internal density of an eigenstate of the internal Hamilto- 
nian Ti nt + U with eigenvalue E int . We have 



IV. EXAMPLE: TRANSLATIONAL MOTION 

An illustrative example is useful at this point. 
Consider the case, already mentioned above, of a 
translationally-invariant Hamiltonian, where we use the 
components of the CoM coordinate vector R as collective 
coordinates: Qi(r) = x/N, Q 2 (f) = y/N, Q 3 (r) = z/N. 
As mentioned above, translational symmetry allows to 
write the energy as a functional of the internal density 
without further consideration for collective motion. How- 
ever, applying the formalism of the previous section to 
this case is useful, since all quantities entering the collec- 
tive Hamiltonian can be derived analytically. 

Using these definitions, 



\f(R)\ 2 = (P(R)) = - S in 2 (K.R), (73) 
with the phase 

9(R) = { ° fol ° - ^ ' ^ ~ 2m7T < n (74) 
{ ' \ 7T for vr < K ■ R - 2mir < 2ir ' 

for integer m. The wave function minimizing the en- 
ergy for this choice of / has ^ cm (R) = f(R). Note that 
J M , Eq. (|57p is, in this case, the average momentum of 
the system. This quantity is zero for this state, thus 

MR)=0. (75) 
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We can derive the collective potential by using the 
above and Eqs. ([73]) and ([74j) in Eq. ([60]): after some 
trigonometry and much cancellation, 



the 4>i(q,r) being a set of orthogonal single-particle 
orbitals. The notation {4>i(q)} p q restricts the varia- 
tional domain to sets of orbitals satisfying 



U(q) 



TT 2 K 2 

2N 
~2N 4 



^26 2 (K ■ R-mir) 



N 



E. 



(76) 



Here, 8 2 {x) refers to the pseudo-distribution which 
yields zero for functions having a node at x = 0, and 
infinity otherwise (which is the result obtained by mak- 
ing 9(R) vary smoothly from to it on an interval whose 
width is then taken to zero) . The expectation value of U 
is thus 

J d n q f* (g) U{q)f{g) = E-^ = E int . (77) 

In the Schrodiiiger equation, this term will constrain 
the cwf to have nodes at K ■ R = mir, i.e. the same as 
our original choice for /. Since the phase 9(R) which 
introduces this term is indissociable from the choice of 
/ (the only reasonable phase choices are ones that make 
/ continuous), the expectation value of this operator, in 
fact, vanishes for any trial cwf. It can thus be dropped 
from the collective Schrodingcr equation. Similarly, V M 
vanishes except for a similar singularity at the nodes of 
the wave function, and the same observation applies. 

Finally, U Ta is in this case proportional to and 
thus negligible. The collective Schrodingcr equation we 
obtain thus involves a kinetic term with the mass of the 
nucleus and a constant potential equal to the internal 
energy. The cwf / chosen initially is a trivial solution 
with eigenvalue and expectation value equal to the total 
energy of the system. 



INTRODUCING ORBITALS: 
SCHEME 



KOHN-SHAM 



In this section I attempt to introduce single-particle 
orbitals in the formulation of the previous sections. This 
introduction is intended to be "minimal" , i.e. as simple 
as possible — other formulations could be envisioned. 

Let us start by writing the collective potential of 
Eq. (JH3I) as 



■W M [p q ]+U ic [f,d}(q), (78) 



U[f,d](q) = T s [p g ] 



where, for convenience, we define p q (r) = d(q,f) and 
T s is the usual Kohn-Sham kinetic energy functional 



(80) 



Note that we use the unmodified mass of the parti- 
cle in the kinetic operator. It is common, in nuclear 
functionals, to include a CoM motion correction, either 
by simply multiplying the particle mass by (1 — 1/N) 
or also including, in addition, the two-body part of the 
internal kinetic-energy operator in the energy. Bear in 
mind that the decomposition ([75]) and (|T9"|) is merely a 
choice, which should be judged on its practical merits. 
From first-principle arguments, the use of such a CoM 
correction in a functional of the internal density is not 
required 0, [5(| ■ Note that we built our formalism start- 
ing from a many-body ground-state wave function with 
a vanishing CoM kinetic energy, and do not need to re- 
move the latter. The Slater determinant formed by the 
KS orbitals has a non-vanishing CoM energy, but it is 
only a theoretical auxiliary. Moreover, since its explicit 
particle-number dependence breaks size consistency, such 
a correction is undesirable in applications to reactions 
[5l| or fission [12,ll3l> hence its omission here. 

The term U ayA captures the contribution from the ex- 
ternal potential, 



U 



' cxt [ Pq \^ I d 3 rv cxt (r)p q _(r), (81) 



while U lc is the interaction and correlation contribution 
not accounted for by the previous terms. In canonical 
electronic DFT, interaction and correlation terms in the 
functional are further split into Hartree and exchange- 
correlation terms. We shall keep the formalism more 
compact and general with respect to the form of the in- 
teraction by omitting this step. 

The density d(q, r) can be optimized by minimizing the 
energy with respect to the orbitals with a normalization 
constraint, as well as a constraint on the average value 
of the collective coordinates to satisfy Eq. (|24|) . 



S E - e k {q)(qk\qk) -Y^J^QM 



dE 



which, per Eqs. <j63j), ([78]) and ([79]), yields 



0, (82) 
(83) 



T s [pq\ = ^ mm 



N 

~2j d'rJ^^k^AM^r) 



(79) 



- + v ext (r) + v s (q;f) 



A M Q M (rO cf> i (q;r) = 0. (84) 
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where we have redefined the singlc-partice energy as 
6k(q) = |/(<z)|~ 2 efc(g) and the Legendre multiplier as 
A^ = 1/(9)1-"% 

The auxiliary potential v s (q,f) is 



v s (q;r) = \ f(q) 



- 2 6E[f,d] 



8d(q, r) 
\.f(q)\- 2 I d n q' f*(q') 



(85) 
(86) 



l„,5A^{q') 5W c {q') 
-o'—r-n — =^d, 1 



2 " 8d{q,f) 



Sd(q, r) 



/(?') 



This formulation raises the usual problem of non- 
interacting v-representability [44j . i.e. of the existence 
of a map between v(q,r) and d(q,r) subject to Eq. ([53]) 
at each q, which presents itself in the same way in the 
present formalism as we use the KS kinetic-energy func- 
tional T s . 

Per the Hohcnberg-Kohn Theorem applied to the non- 
interacting system, Eqs. (|75|) and ([M]) unambiguously 
define a unique potential V s (q,r) and a unique set of or- 
bitals and associated energies which, for each value of q, 
are functionals of p q {r). These orbitals are labelled by an 
index k which can be smaller or greater than N. In the 
following we use i as an index on the first N (occupied) 
orbitals, a as an index on unoccupied (or virtual) orbitals, 
while k indexes the whole basis. This in turn allows to 
introduce orbital-dependent terms in the functional. Let 
us use this possibility by expressing the collective mass 
as 



-V [/, d] (g) = A 1 ^ [d] (q) + A% [/, d] (q) . (87) 

In the last expression, A]^ , is the Inglis cranking- 
formula collective mass [13, Hj] , while A^ v is the remain- 



ing interaction-correlation component. 



A ln (g) 

where 



[B in (g)}-\ 

„ (£i|<°M(Ht4 \gi) fon . 

% " {e a {q)--eMf~ ' ^ 



dv{q, r) 
dq» 



(90) 



(qk\vM\ql)= d 3 rtt( qi r)v^\r)Mq,r). (91) 



Orbital-dependent terms have been extensively used in 
quantum-chemistry applications of DFT, first as a way of 
replacing non-local exchange terms [HI, [5(| , see also [j57| , 
then as a tool for introducing explicit correlations in the 
functional through perturbation theory [58l . [59j . Let me 



refer to the reviews in Refs. [8j, [60J, [61| and simply use a 
straightforward generalization of their main result (make 
all quantities depend on q and add a factor |/(g)|~ 2 to 
the definition of v s ) to derive the contribution vi n (<7, r) of 
A l ^ v to the auxiliary potential. This can be obtained by 
inverting the optimized effective potential (OEP) equa- 
tion, 



l/(?)r / df Xs(q;f,nvin(q)n = A In (< ? ;r), (92) 



where 



Ai„(g;r-) = ^|-| 



df 



0l(q;r)Gk(q;r,r') 7-^ 



\Mq;?)f 



5E ln 



(93) 



8e k (q) 

and E\ n is the cranking contribution to the energy, 

Em[f,d] = -\( d n qf*{q)Y,^ (A%[d](q)d„f(q)) , 

(94) 

Xs being the static KS response function and G k the 
KS Green's function, 



3 ( 2 ; r, 7) = - ]T 4>l (9; r) G k ( £ ; r, 7) <f> h ( £ ; f) 



c.c. 



^ ei(q)-e k (q) 



(95) 
(96) 



l=tk 



The functional derivatives of the orbital-dependent en- 
ergy with respect to and e k , using t as a shorthand 
for either, are 



-ST=2j d n q_r{q)d[{B^)-\q) 
x{B l *)-\q)dj{q) 



(97) 



where the derivatives of B are given by 



' ^4^(r)^(g;f) 



(fc < AT), 



2 ^ 



E 



(gq|«£^jgfc) 
* (s a (q)-e k (q)r 

( £fc (g)- £l (g)) 3 

x vMtfMq;?) (k>N), 



(98) 
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and 



Sek(q) 



= 6 



E(gk\Qn\ga)(qa\Qv\qk) 

E(qi\Qn\qk)(qk\Qv\qi) 
, (e k (q)-eM) 4 



(k < N), 

(99) 
(k > N), 



This completes the sets of equations needed to solve 
for energy-minimizing f(q) and d(q,r). The solution 
should proceed by alternating between Eqs. ([84")) a t a 
set of points in g-space and |65|) . starting from an initial 
guess and iterating to convergence, using standard tools 
for the self-consistent solution of Kohn-Sham equations. 
This is obviously more involved that standard DFT, ow- 
ing to the multiplication of the computational load by the 
number of q mesh points, but this problem appears rela- 
tively easy to treat with parallel processing, as only the 
cwf and fields have to be communicated between neigh- 
boring points, as well as densities unless one assumes a 
lack of dependence of U{q) on densities at q' ^ q. The 
largest data sets, i.e. orbitals, stay local. 

In the treatment of Hill- Wheeler or collective Hamil- 
tonian equations, feedback from collective motion to 
the single-particle "mean field" is usually ignored [281 ]. 
Above, we have a recipe for going beyond that approxi- 
mation, which would be of interest for the description of 
rotational bands in collective nuclei [6^, [HI, as well as 
the dynamics of fission processes ■ 

If we nevertheless neglect the feedback from collective 
motion, the formalism can be put into a more conven- 
tional form. Assuming the .A-term from the equation for 
v s (q,r), Eq. (|87|) to be negligible, the latter then only 
contains a functional derivative of U. The collective po- 
tential U(q) generally depends on d(q' , r) and f(q') for all 
q', due to the coupling between different slices ^(q) in the 
many-body Schrodinger equation. If we further assume 
that this coupling is weak, we can reduce the dependence 
of U on the density at q, i.e. 



U(q)[f,d]=U(q)[p q ] = U[ Pq 



(100) 



T s [p q }+U ext [p q }+U ic [ Pq }. (101) 



where the second equality in Eq. (|100p is justified by 
the fact that the value of q is encoded in p q (f) = d(q, r), 
Eq. pljl. We then obtain 



v s{q,r) = — — . 

OPq 



(102) 



Then, Ui C [p q ] plays the role of the usual Skyrme, Gogny 
or relativistic functional, which can be used in the stan- 
dard way, with an independent, a posteriori solution of 
the collective Schrodinger equation. This form, however, 
relies on the assumptions above. 



More generally, a practical application of the present 
theory requires a parameterization of U lc and A™ v . Let 
me stress the latter: the cranking formula used above is 
by no means assumed correct by itself. It could be re- 
placed by a term derived usin g ad iabatic time-dependent 
Hartree-Fock (ATDHF) [H ED, M, HI or the Gaus- 
sian overlap approximation to the GCM (GCM-GOA) 
J2(| UK- 40]. Even then, such terms rely entirely on the 
single-particle orbitals at each position in the collective 
space and cannot be expected to correctly reproduce the 
physics of the underlying correlated many-body state. 
The collective mass A^ v is an integral part of the func- 
tional, and, in the absence of a rigorous ab-initio deriva- 
tion, it is perfectly reasonable to parameterize it and ad- 
just the parameters to experimental data. 



VI. SUMMARY AND OUTLOOK 

The energy of a many-body system is expressible as 
a functional of a generalized density, which extends the 
concept of the local particle density to include a depen- 
dence on coordinates describing collective motion of the 
particles. The generalized density can be decomposed 
as the product of the square of a cwf and a density pa- 
rameterized by the collective coordinates, which is al- 
lowed to break spatial symmetries of the Hamiltonian. 
By decomposing the kinetic contribution to the energy, 
the functional can be written into a form that allows to 
write a Schrodinger equation for the cwf. When the col- 
lective coordinates are chosen to be the components of 
the inertia tensor of the system, the collective Hamilto- 
nian takes the form of a generalized Bohr Hamiltonian. 
Single-particle quantum effects can be reintroduced with 
s.p. orbitals determined from a s.p. potential deduced 
from the parameterizations of the collective mass and po- 
tential. With the assumption of weak coupling between 
single-particle degrees of freedom at different points in 
the collective space, the functional can be reduced to 
a form similar to current nuclear Energy Density Func- 
tional augmented by a collective Hamiltonian; equations 
have been derived for going beyond this scheme and op- 
timizing s.p. and collective degrees of freedom simulta- 
neously. 

Extending the formalism to supernuid systems 
presents a significant challenge, and it is likely that a sim- 
pler scheme can be found for this particular case. Deriva- 
tions have been carried out ignoring spin and isospin 
degrees of freedom; reintroducing these [ 50f . as well as 
introducing spin and kinetic densities [671 |68[ ] , seems to 
pose no major obstacle. Developing a time-dependent 
version of the theory presented here in the vein of Refs. 
jr39l |70T ] could prove useful for the treatment of excita- 
tions of deformed and highly collective nuclei, as well as 
nuclear reactions [7l| . Finally, it would be interesting to 
derive from first principles, for a few nuclei, the collective 
Hamiltonian proposed here. This is made difficult by the 
TV-body operators involved, but probably feasible using 
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a many-body method that uses the 3V-dimcnsional co- 
ordinate representation nativel y, s uch as Variational or 
Green's Function Monte-Carlo 1721173^ 
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